Polarization operator for plane-wave background fields 
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We derive an alternative representation of the leading-order contribution to the polarization op- 
erator in strong-field QED with a plane-wave electromagnetic background field, which is manifestly 
symmetric with respect to the external photon momenta. Our derivation is based on a direct evalua- 
tion of the corresponding Feynman diagram, using the Volkov-representation of the dressed fermion 
propagator. Furthermore, the validity of the Ward-Takahashi identity is shown for general loop 
diagrams in an external plane-wave background field. 
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I. INTRODUCTION 

The most precise calculations known so far in physics 
are provided by QED. The reason for this is the smallness 
of the fine-structure constant a — e 2 / ' (Aire^hc) sw 1/137, 
which allows to use perturbative techniques at sufficiently 
low energies (e is the electron charge). The most promi- 
nent example is probably the electron (/-factor, for which 
experimental and theoretical results have been matched 
on the record accuracy level of parts per billion [l|. 
To achieve this outstanding precision, the corresponding 
theoretical calculation included Feynman diagrams with 
up to four loops. 

A quite different situation is encountered for QED 
with electromagnetic background fields. A source for 
strong electromagnetic fields are modern laser systems. If 
spatial focusing effects are sufficiently small, laser fields 
can be well approximated by plane-wave fields. For a 
plane- wave field we obtain, besides the fine-structure con- 
stant, a second gauge and Lorentz invariant parameter 
£o = \e\Eo I '(mcoJo) , where Eq and luq are the peak electric 
field strength and central frequency of the plane wave, 
respectively (m is the electron mass). If Co ^ 1 the in- 
teraction between electron and positrons with the laser 
field must be taken into account exactly. For optical 
lasers (photon energy fkuo sa 1 eV) this happens already 
at intensities of the order of 10 18w /cm 2 . More precisely, 
we can generally still treat the interaction of the elec- 
trons and the positrons with the quantized radiation field 
perturbatively as in vacuum QED (QED without back- 
ground fields) , but must include the dependence on £o to 
all orders if this intensity is exceeded. 

Another important scale is the so-called critical field 
E cr = m 2 c 3 /(h\e\) = 1.3 x 10 16 V/cm, which corre- 
sponds to a peak laser intensity of I cr = eoC-E 2 r = 
4.6 x 10 29w /cm 2 . A constant and uniform electric field of 
this strength can, in principle, produce electron-positron 
pairs from the vacuum [2j-|4|. 

The current laser intensity record (in the optical 
regime) is given by 2 x 10 22 w /cm 2 !5| and future facilities 



envisage even intensities of the order of 10 24 — 10 25 w /cm 2 
[fJ-Q. Thus, the non-perturbative regime (in £ ) can be 
entered with presently available laser systems, and even 
the critical field can be reached in the rest frame of an 
ultra-relativistic particle (e.g. a ~ 1 GeV electron j9[). 

So far only one experiment has been carried out to 
probe strong- field QED effects using laser fields piJ, [HI . 
However, this is expected to change in the near fu- 
ture and, correspondingly, the experimental progress has 
stimulated many theoretical investigations during the 
last years [13 - [46( | . For a more detailed overview the 
reader is referred to the review [47j |. 

In contrary to vacuum QED, calculations with a plane- 
wave background field have not been carried out beyond 
the one-loop order (for constant-crossed fields higher- 
order calculations have been performed, see e.g. |48l - t50| V 
This can be attributed to the fact that already diagrams 
with just a few propagators correspond to quite compli- 
cated expressions. It is therefore of general interest to 
investigate new techniques, which have the potential to 
make also the calculation of complicated diagrams trace- 
able. 
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FIG. 1. The Feynman diagram corresponding to the leading- 
order contribution to the polarization operator 7"*" (91,92) in 
a plane-wave background field. The double lines represent the 
Volkov propagators for the fermions, which take the external 
field into account exactly [see Eq. (|19p ]. The vertical dashed 
line links the polarization operator to the pair-production di- 
agram due to the unitarity of the S'-matrix. 



Here we present a new derivation of the first-order 
contribution to the polarization operator given in Fig. [T] 
|5ll I52J . In [5l[ an operator approach similar to the one 
introduced by Schwinger [J] has been used. We show 
here how the diagram in Fig. [1] can be evaluated directly 
using the Volkov-representation of the dressed propaga- 
tors. This approach has the appealing feature that it is 



very similar to the calculation techniques used in vac- 
uum QED. Our final result (which is equivalent to the 
one in [51|) has the interesting property that it is man- 
ifestly symmetric with respect to the external photon 
four- momenta q\ and qi- Furthermore, we prove the va- 
lidity of the Ward-Takahashi identity |53|, [54[ for general 
loop diagrams in a plane-wave background field. The 
calculation techniques employed here are expected to be 
useful also for other higher-order calculations with strong 
plane- wave background fields. 

The polarization operator itself is of central impor- 
tance, since it determines the properties of a photon in- 
side the background field via the Schwinger-Dyson equa- 
tion for the exact photon propagator [5l|, [55[ . As a con- 
sequence, a plane-wave field acts as an active medium, 
e.g. the photon obtains a mass and has a non-trivial dis- 
persion relation. Furthermore, due to the unitarity of the 
S- matrix, the imaginary part of the polarization operator 
is related to the total photo-production probability for an 
electron-positron pair [5a, [57[ . The polarization operator 
is also required for the calculation of radiative corrections 
to elementary processes like non-linear Compton scatter- 
ing or pair production. The significance of the photon 
polarization tensor can also be inferred from the ongo- 
ing effort to calculate it for different field configurations 

The paper is divided into three parts. In section IPTl the 
general framework of strong-field QED with plane-wave 
background fields is described. The actual calculation is 
then presented in section IIIII Finally, we show how our 
results are related to those obtained in [51( and discuss 
the special cases of a constant-crossed field, a linearly 
polarized field in the quasi-classical approximation and a 
circularly polarized, monochromatic field in section IIV1 



A. Vacuum QED 

To obtain a quantum theory of electrons, positrons 
and photons, both the spinor if) and the photon field A 11 
are promoted to operators with canonical (anti-) com- 
mutation relations (alternatively, the functional integral 
formalism can be employed). Once derived in either 
way, the S'-matrix element for a given process can be 
calculated perturbatively using Feynman rules. In vac- 
uum QED the starting point for the perturbative expan- 
sion are the solutions of the free Dirac and wave equa- 
tions [Eq. © with A* = and Eq. © with J» = 0], 
An electron with a given four-momentum p 1 ' = (e,p) 



(e > 0,p 2 = 
wave solutions 



% 



can then be described by the planc- 



6 Upi 



m)u p = (4) 



and the corresponding propagator is given by 



iG{x,y) = i 



d 4 P 



-ip(x-y) _ 



i> + in 



p 2 — m 2 + iO 



(5) 



For a photon with polarization four-vector e M and four- 
momentum fc M = (u>, k) (w > 0, k — 0) we obtain the 
following wave-function 



f 
2w 



K = -4=e- tkx e^ 



(0) 



and in the Feynman gauge the photon propagator is given 
by 

-iD^x-y) = -i [^e-^'^J^-. (7) 



(2tt)< 



k 2 +i0 



II. STRONG-FIELD QED 

QED is described by the following Lagrangian density 

MM 



c 



QED 



4>(i0-m)il>- \T^ V - mhy&A*, (1) 



where and A^ are the Dirac spinor field and the electro- 
magnetic four-vector potential, respectively, and T^ v = 
d^A v — d v A^ is the electromagnetic field tensor [from 
now on we will use natural units h = c = 1 and Heavisidc- 
Lorentz units for charge, a = e 2 /(47r), see appendix lAl 
for more details] . The equation of motion for the spinor 
field ip, which follows from Eq. (JXJ) , is the Dirac equation 



(ip ~ eA — m) -0 = 0. 



(2) 



Correspondingly, we obtain for the photon field A^ in 
Lorentz gauge (d^A^ = 0) the wave equation 



dPdpA" = eJ», 



J^ = tfj^tfj. 



(3) 



Finally, the interaction between electrons, positrons and 
photons is represented by the interaction vertex 



-le 



d 4 x ■■■ 7 P 



(8) 



where the dots indicate that the vertex is always con- 
tracted with two fermion and one photon line. Thus, 
one can move the exponential functions from the ex- 
ternal lines and propagators to the vertex. After the 
space-time integrals associated with the vertex arc taken, 
momentum-conserving delta functions are obtained. For 
a more detailed discussion see e.g. |55ll60j. 



B. QED with background fields 

Vacuum QED has been tested to a very high precision, 
because, due to the smallness of the fine-structure con- 
stant, a perturbative treatment is adequate in most situa- 
tions. However, for very strong external electromagnetic 
fields the (conventional) perturbation series breaks down. 



Modern laser facilities provide a source of such strong ex- 
ternal electromagnetic fields. A laser field is described by 
a coherent state of the photon field that effectively leads 
to the splitting A^ = A^ ad +A^ in the Lagrangian. While 
A^ ad represents the quantized radiation field, the classi- 
cal field Af 1 describes the laser field [61l - l63l ]. Working in 
the so-called Furry picture [6J| the classical background 
field A^ is taken into account exactly and only the radi- 
ation field is treated perturbatively. The starting point 
for the perturbative expansion of the S'-matrix is then 
the solution of the interacting Dirac equation @ with 
the replacement A^ L — > A 1 - 1 . Since photons have no self- 
interactions at tree level, the photon wave functions and 
propagators are left unchanged by the background field. 
A more detailed discussion of strong-field QED can be 
found in [47], |5CJ, |6l|, |65| and in the references therein. 



C. Plane- wave fields 

In this paper we will consider only plane- wave external 
fields A^ . This means that the field tensor F^ v depends 
only on the plane-wave phase <j) = kx, where /c M is a mo- 
mentum four- vector which characterizes the plane wave. 
Such a field can be described (in the Lorentz gauge) by 
the following four- vector potential 0, [5l|, |66| 



A»(kx) = a1ipi(kx) + a^ 2 {kx) 



(9) 



where k 2 — ka,i — a\a 2 — and where tpi(kx) are 
arbitrary functions, restricted only by the physical re- 
quirement that the external field is of finite extend (e.g. 
ipi(±oo) — ip'^zkoo) — and we also assume that the 
field vanishes fast enough at infinity). Furthermore, we 
adopt (without restriction) the normalization condition 
\ipi(kx)\ , |V>i(fcaO| ^ 1 which means that the strength of 
the field can be characterized by the Lorentz invariant 
parameters 



6 



f 



(10) 



which we will call the classical intensity parameters. It 
turns out that the plane-wave field must be taken into 
account exactly if & > 1 [47]. Modern laser facilities 
can easily reach this non-perturbative domain, e.g. in [5[ 
£i ~ 100 was obtained. 

The field tensor corresponding to the four-potential in 
Eq. (J9J is given by 



F^(kx) = /rmkx) + frmkx), 



where 



f„P 



ir = fc'x - ax 



(ii) 

(12) 
(13) 



It is also convenient to introduce the integrated field- 
strength tensor 



r u (k X ) 



#'i^"(A 



(14) 



which can be written as 

^(kx) = k lx A u {kx) - k^A^ikx) 
= frMkx)+frMkx) 



(15) 



in the Lorentz gauge [we will use 3*f" = ^ u (kx) inter- 
changeably to denote the argument]. Both F^ v and y" 
have the important algebraic property that successive 
contractions of more than two tensors vanish and their 
square is proportional to fc^fc", e.g. 









i=l,2 



a 2 ipi(kx)i>i(ky). (16) 



If the background field is a plane- wave field, the Dirac 
equation [Eq. © with A^ — > A**] can be solved analyt- 
ically [67J. The corresponding so-called Volkov-solution 
with the boundary condition \P P — > ip v if kx — > — oo [see 
Eq. flU] can be written as [il. l55l I68J 



*„ 



£ J p,x'U'pi 



-t^p.x — 



1 + 



,;) 



2kp 



JS p {x) 



(17) 



where the phase is given by the classical action 



S p (x) 



-px - 



k:r 



epA(4>') 
kp 



2 A 2 



2kp 



#'. (18) 



Despite being essentially semi-classical, the Volkov solu- 
tions are exact solutions of the interacting Dirac equa- 
tion. Correspondingly, the dressed propagator (which is 
the Green's function of the interacting Dirac-equation) is 
given by 



iG(x,y) = i 
where 



d A p 

(2^ Cp 'V-m 2 



En 



i0 



:E V 



^p 7 x — 



zA-(kx)$ 



2pk 



-iS p {x) 



(19) 



(20) 



Thus, in comparison with the vacuum case, the plane 
waves are replaced by the Ritus .Ep-functions, which de- 
pend non-trivially on the plane- wave phase kx. However, 
they also form an orthogonal and complete set |48J 



d 4 p 
(2tt)* 



&p,x-&p.x / ' [& $ )•) 



(21) 



d 4 xE p , iX E^x^(2Ti) A 5 A {p' -p). 



The _E p -functions convert the dressed momentum into the 
free momentum 148] 



[i$ x - eA{kx)]E^ x = E PtX f, 



-idjfEp^jp - eE p ^A{kx) = 



(22) 



J p,x 



(these identities hold only if the derivative acts solely on 
E PtX and E p>x , respectively). 



D. Dressed vertex 

To obtain Fcynman rules in momentum space, we can 
proceed analogously as in the vacuum case and move the 
-Ep-functions to the vertex [68|]. Correspondingly, we de- 
fine the dressed vertex by 

V(p',q,p) = -«■ / ' d A xe-^E p , jX YE p ^. (23) 



Working in momentum-space, the only difference be- 
tween vacuum QED and strong-field QED is the vertex 
we have to use [i.e. the free vertex in Eq. ([8]) is replaced 
by the dressed vertex in Eq. (|2"3"j) ]. Using the relations 
given in appendix [C] we can write the dressed vertex as 



T p (p',q,p) = -le I d 4 x [^G pp (kp',kp;kx) 

+ ijnl 5 G pp {kp',kp;kx)]e zSr{p '' q ' p - x \ (24) 
where the phase and the coupling tensors are given by 



Sr(p',q,p;x) = -S p >(x) -qx + S p (x) 

r kx . r«n..r>'3/"'l,/l 
(p -q-p)x- 



ep P p'JB p 



e 2 (kp— kp') 



(kp)(kp') 



TSPM^m , (25) 



2(kp) 2 (kp') 



G pp (kp', kp; kx) = gw + Gx$» + Gtf%*, ( _ 

*w ( 26 ) 



G% p (kp', kp; kx) = G 3 $* x 



Gi = 



kp + kp' 
-e — — ^-r, G 2 = 



G 



2kp kp' 2kp kp' ' 

kp — kp' 



(27) 



2kp kp' 



(note that G\ and G2 are even in the permutation kp -H- 
kp' while G3 is odd). We point out that the expression 
given in Eq. (|24|) is manifestly gauge-invariant, since it 
depends on the external field only through the tensor $ pi/ 
[68| . Furthermore, if we subtract the vacuum vertex and 
consider the quantity 



T p (p',q,p)+i el p 



d*xe i( - p '- q - p \ 



(28) 



all integrals are properly convergent. 

In position space the dressed propagator in Eq. (|19l) 
can be interpreted such that the electron (or positron) 
continuously interacts with the external field during its 
propagation. Examined in momentum space, we can also 
visualize the influence of the external field as a modifica- 
tion of the coupling between the photons of the radiation 
field and the charged particles. From Eq. (|2"4"|) we see that 
beside the modification of the photon vector current in- 
teraction we also obtain a coupling to the axial-vector 



current inside the plane- wave background. This is pos- 
sible, since the external field provides the pseudo-tensor 

Since the external field depends only on the plane- wave 
phase cj> = kx, it is useful to use light-cone coordinates, 
which are discussed in appendix [Bl We can then always 
take the integrals in dx + and dx^- in Eq. ([24|) and obtain 
momentum conserving delta functions in three of four 
light-cone components 



Sl-^tf-p-q), 



where we used the notation 



S { -'^(a) ^ 5(a-)S(a 1 )5(a 11 ) 



(29) 



(30) 



for a general four- vector a p . Thus, the four- momentum 
is only conserved up to a multiple of the plane- wave four- 
momentum k p at each vertex. 



E. Ward-Takahashi identity 

The Ward-Takahashi identity [53|, [5J] is a direct conse- 
quence of the gauge-invariance of QED, which becomes 
Mticular transparent in the functional integral approach 
l69| . Diagrammatically, it is a functional relation 
for Feynman diagrams (in momentum space), where the 
polarization four-vector of an external photon leg is 
replaced by the corresponding momentum four-vector. 
In [70] a perturbative proof of the Ward-Takahashi iden- 
tity in vacuum QED is given. We will show now how this 
proof can be extended to electron-positron loops inside a 
plane-wave background field. 




FIG. 2. Closed electron loop with n dressed vertices and 
propagators. 

The starting point is the following algebraic identity 
for the dressed vertex 16811 



q p TP(p',q,p) 



where 



i' - m)I(p', q,p) - I(p',q,p)(p - m), (31) 



I(p',q,p) = -ie / d 4 xe- iqx E p ,, x E p , x . (32) 



To verify Eq. ([3T]), we use E q. ((22J) and note that the 
identity 



d i xid lt [E p >0'fe- i9 *E p , x -e i W- 9 - p ')''f] =0 (33) 



holds (since the external field is supposed to vanish at 
infinity). 

Typically, (p 1 — m) and (p — m) in Eq. (|3ip cancel 
an adjacent propagator and the associated momentum- 
integral can be taken using the relations 



4„// 



I^,q , ,p")r^(p",q,p , ) = -ieT^(p,q + q , ,p'), 

T»(p,q,p")I(p" ,q' ,p') = -teT"(p,g + g',p , ) > 

(34) 



(2^) 4 



which follow from Eq. (plj). Using Eqs. (|3Tj) and (|34|) 
we can simplify diagrams which contain dressed vertices 
contracted with the corresponding photon four-momenta. 
As an example, we consider now a closed electron loop 
which contains n dressed vertices and electron propaga- 
tors (see Fig. [5]) . The i th propagator of such a loop to- 
gether with its adjacent vertices is given by 

• • • r«( Pi _i jft ,p i )— ?— r^Cpi, q i+ up i+ i) ■ ■ • (35) 



(the electron four- momenta pi are integrated out). If 
we insert now a vertex (contracted with its photon four- 
momentum) at this propagator, we obtain 

■■•T^ipi-uq^pi)- q^(pi,q,p')x 



r^(p',q i+1 ,p i+ i)--- (36) 



and, by using Eqs. (|3"Tj) and ((Ml) , we find that this is 
equivalent to 

■■■T^(jpi-x,qi + q,pi)- T^ i+1 (pi,q i+ i,pi + i)--- 



T^ipi-i^upi)- r w+1 (pi,g i+1 +q,p i+ i)--- 

p. — m 

(37) 

(in the first line we have changed the name of the inte- 
gration variable from pt to pi). Thus, the insertion splits 
the diagram into the sum of twice the original diagram 
with the additional four-momentum q added once at the 
adjacent vertex before and after the insertion. If we sum 
now over all possible insertion points of the loop, we ob- 
tain zero, since all contributions cancel (as in the vacuum 
case |70|). 

We point out that the above discussion is shortened, 
since possible issues arising due to the renormalization of 
the theory were not addressed (in general, the validity of 



the Ward-Takahashi identity may be spoiled by anoma- 
lies [71j). In this paper, however, we are mainly inter- 
ested in modifications induced by the background-field, 
which turn out to be finite. Thus, subtleties arising from 
manipulations of divergent integrals can be addressed as 
in vacuum QED. 



III. POLARIZATION OPERATOR 

A. General expression 

The leading-order contribution to the polarization op- 
erator V J ' l '(qi,q2) for plane- wave background fields (see 
[55| . §104) is determined by the diagram in Fig. [TJ This 
diagram corresponds to the following expression 

T^( qi ,q 2 ) = J*^ftrr»(p\q u p)x 

(p + m) 



p 2 — m? + iO 



r»(p,-q 2 ,p') 



pll _ m 2 _|_ jQ 



(38) 



and T» v = iV» v (see j55|], §115; [72|)- We note that 
T^ v (qi,q 2 ) is divergent, but if we write 

T^( qi ,q 2 ) = [T^( qi ,q 2 )-T^ ( qi ,q 2 )] 

+ %(?!,<&) (39) 

the first part is finite [5l| and the regularization of the 
vacuum contribution is well known [55l |60( • In the fol- 
lowing, we will focus on the first part which contains only 
the corrections induced by the external background field. 

To determine the expression in Eq. ([38]) we have to 
insert the dressed vertex given in Eq. (|24[) [we will de- 
note the vertex integrals associated with r M (p',<7i,p) 
and T v (jp, —q 2 ,p') by d A x and d 4 y, respectively]. We then 
obtain for T» v (qi,q 2 ) 



T^(q x ,q 2 )=±{-ief j ' ^f J ' d*xd* 

itr[-r 



yx 



(p 2 — m 2 + i0)(p' 2 — m 2 + iO) 



JS T 



(40) 



(the prefactor 1/4 is included explicitly for later conve- 
nience), where the phase reads [see Eq. (|24l) ] 



iS T = i(p' -p- qi)x + i(p - p' + q 2 )y 



kx 

i I d<t>' 

Iky 



ePv,p' v ^ v , e 2 {k P -kp') 



/\ 2 P/J,P v -0 



(41) 



(kp)(kp') 2{kp) 2 (kp') 2 ' 
and \ tr [ • • • 1 in Eq. (|40p can be calculated using the 



relations given in appendix [C] 

i tr [ 7a a°" + i lal b b a »] (f + m)x 

x [ l0 c^ + i lpl b d^} (/ + m) 
= m 2 [(a<*"0 + (6°**0] 
+ (wO^a") - (pp'X^O 
+ (p Q a^)(^c^) + (p' a a^)( W c^) 

- ( Pa b a n(p'^n - (p'^mppd^) 

- e paa0 p p p^(a a ^d^ + b a ^c^), (42) 



where 



a a M = Q.*^ fcp . fca ,) ; c ^ = G^(fep, fcp'; %), 
feQM = G5 M (/cp', ftp; to), d' 3 " = Gf " (ftp, V; *y). 



B. Evaluation of the integrals 

Working in light-cone coordinates (see appendix [B]) we 
can take all space-time integrals except of those in dx~ 
and dy~ , and obtain the momentum-conserving delta 
functions 

(2w) 6 S^^ (p'-p- q) 5(-^ ( qi - q 2 ). (44) 

Here and in the following we write q^ if q p and q% can 
be used interchangeably due to the above delta function. 
Successively, we can take the integrals in dp'~ and dp ,J - 
(for simplicity we will continue writing p' and identify 
p' = p + q for the components — , J_). 

It is now convenient to introduce the two four-vectors 



A? = 



/i 9. 



kqyj-a\ 



, A 



fi _ J2 



£1X11 



kqv~ 



(45) 



which obey A^Aj = — <5y, kAi = qiAj = and 

/rAi„ = -?*"&, / 2 ^A 2 , = -f ft^ 2 . (46) 

They allow us to write the remaining phase as 

iS T = i{p' -p- qi) + x~ 

+ i{p-p' + q 2 ) + y~ +ip\ + iA, (47) 



where we defined 



A" 



A= - 



> i=l,2 "'** 



m(kq) 

'(*p)(*pO ,fr 

m 2 (ftg) y 
2(kp)(kp') ^ 2 .,,„ 



(48) 



Due to the appearance of A^ in A M , it is more convenient 
to use modified light-cone coordinates from now on [see 
Eq. (jBlip . the calculation so far is independent of this 



choice] . In modified light-cone coordinates we obtain the 
convenient relations 



-Ll-L 



pA = — p A 



q X =0, 



P' X =P X , 



(49) 



which simplify the algebra considerably. 

If the preexponent would not depend on p + and p' + , 
both integrals could now be taken. We therefore intro- 
duce the pro per-time representation of the scalar propa- 
gators [J, [63 



1 



1 



m 2 + zO p' 2 — m? + iO 



dsdt x 



x exp [i(p 2 - m 2 + iO)s + i{p' 2 - m 2 + iO)t] . (50) 



(43) In the following, we will drop the pole-prescriptions «0 



and keep the replacement m 



— iO in mind. Fur- 



thermore, we add the source terms ip^j^ + ip'nj'* 1 to the 
phase, which allows us to make the replacement 



t — * (-0 



f 



H) 



(51) 



in the trace. Now, the preexponent depends only on p~ 
(through ftp and ftp'). Taking the derivatives with re- 
spect to the sources out of the integrals, we can take the 
integrals in dp + , dp' + which results in the delta functions 

(2tt)%- - x- - ^- t (2stq- - tj~ + aj'-)] x 

x (2Tr)S[2p-(s + t) + 2q-t+j- +j'-}. (52) 

Successively, these delta functions can be used to take 
also the integrals in dy~ and dp~ . To this end we rewrite 
(since s + t > 0) 

S[2p-{s + t) + 2q-t + j- +j'~] 

= 2(h) 5 ^ + 2ih)( 2 ^ t + r+r)] (53) 

(for simplicity we keep writing y~ and p~). In particular, 
we obtain the identities 



kp = —± i [tkq+±(kj + kj')] [ 
kp' =+-L.[ 8 kq-±(kj + kj')] 



(54) 



ky = kx + -jj^{2stkq — tkj + skj'), 
which imply for j = j' = that 
e (s-t)(s + t) 



Gi = 
G 2 = 



2kq st 

e 2 (s + t) 2 



G, = - 



2(kq) 2 st 
e (s + t) 2 



2kq 
where we defined [5l| 

G 

T = S + t, V 



St 



e vt 
2kq /i 

e 2 t 
2(ftg)V 
e t 
2kq \x 



(55) 



s + t' 



st 

S + t 



\r(\~v 2 ) (56) 



[the motivation for these definitions becomes clear in 
Eq. C3)]. 

The remaining part of the phase structure (including 
the part coming from the propagators and the sources) 
is now given by 

iS' T = i[(q+ - q+)x~ + -£- t ql - ^(tcftj - s q 2 f) 

- {p X p X +m 2 )(s + t) 

-(J X +j' X +^ X )p ± +A]. (57) 

Taking the Gaussian integrals in p l and p u we obtain the 
prefactor y-xtl and the final phase is given by 

iS' T = i[(q+ - q+)x- - m 2 (s + t) + ^- t q 2 

-W+tlti+j')*-WWj x2 + A ]> ( 58 ) 
which reads for zero sources (j = j' = 0) 

iS' T = i [(q£ - qt)x~ + fiq 2 

-Tm 2 +Tm 2 Y / &(I?-^)], (59) 



where we defined 
h-- 



1 



2kqn Jky 

ckx 



(MM*'), 



J l = - 



1 



2kq[i J k . 



dtftf (</>') 



fry 



(the prefactor is chosen such that for j = j 
ipi = 1 we obtain Ij = Jj = 1). 

Finally, we can write the tensor T^ v as 

/>oo 

T^(qi,q 2 ) = ~2me 2 8 { --^{q 1 -q 2 ) / ds 

Jo 



(60) 



and 



dt x 



+oo 



dx 



(s+ty 2 4 



tr[ 



y iv e is' T 



j=j'=a 



(61) 



where the expression for i tr [. . .]'"' is given in Eq. (|42|) 
with the replacement in Eq. (|51[) and where the sources 
are set to zero after the derivatives are taken. 

We point out that the two four-momenta q\ and q 2 ap- 
pear asymmetrically in the final expression [see Eq. (159|) ]. 
To remove this asymmetry we shift the x~ integration by 
defining 

z~ = x~ + fxq" '. (62) 

After this shift the phase contains q\q 2 since 

(i 2 - <ii) x ~ + ml = (it - ii) z ~ + Mi 8a- (63) 



Furthermore, we obtain (for j = j' = 0) symmetric rep- 
resentations for the functions in Eq. (|60|) 



(64) 



h = h / dXipi(kz — X/ikq), 



Ji = \ \ dXipf(kz — X/ikq), 



since 



kx = kz — ukq, 

(65) 
ky = kz + [xkq + -A^(skj f — tkj). 



C. Tensor structure 

In principle, the only remaining task is to evaluate 
the two derivatives with respect to j and j' and then 
set j = j' = 0. Despite being elementary, this is the 
most tedious part of the calculation, since the sources 
appear in many places in the final expression. The work 
is considerably reduced if we expand the polarization op- 
erator in a convenient basis [5l| . To this end we note 
that 

qi^ v ( qi ,q 2 )=0, T^( qi ,q 2 )q 2v =0 (66) 



due to the Ward-Takahashi identity (see section Hi E[) . 

Since the four-vectors Aj appear in the phase [see 
Eq. (|47p] and c^Aj = 0, it is natural to introduce the 
two complete sets gi, Qi, Ai, A 2 and g 2 , Q 2 , Ai, A 2 , 
where 



_ k»q* - q^kq _ k^q 2 - qjkq 

«1 — i J «2 



kq 



kq 



(67) 



(fi? 



-9i 



2 Qi 



-g|, Q 4 Aj = 0, g,Qi = 0). Using 



the set including qi for the index /i and the set includ- 
ing q 2 for the index v, seven of 16 coefficients vanish 
due to the Ward-Takahashi identity and we can decom- 
pose T^{qi,q 2 ) as [H| 

T» v = ciA^Ag + c 2 A£A? + c 3 A?A£ 

+ c 7 Q1A^ + c s A1Q^ + c 9 A^ 2 . (68) 

It turns out that also the coefficients ce — eg vanish. If an- 
alyzed perturbatively (with respect to the external field 
coupling) this can be understood from Furry's theorem 
[5l[ [52j . Since a closed fermion loop with an odd number 
of vertices vanishes, only diagrams with an even number 
of external field couplings (eA^) contribute to T^ v . Due 
to gauge-invariance and the fact that T^ v is a tensor, the 
external field can enter only as ^ v (which is linear in 
A^). Since it is not possible to construct a scalar linear 
in 3^ using only the four- vectors q^ , q 2 and fc M , the ten- 
sor structure cannot involve an odd number of the tensor 
3*"" . As a consequence the coefficients c§ — eg (which are 
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linear in A^ and thus in the external field) should vanish. 
We will later see that this is indeed the case. 

The coefficients Cj in Eq. (|68|) can be determined 
by contracting T^ v (91,(72) with appropriate four- vectors. 
Especially, using again the Ward-Takahashi identity, we 
obtain 



ll 



4 



Q lu T^ = ^kuT"", T^Q 2u = ^-T^kv 



kq 



kq 



(69) 



Thus, effectively, we need to determine the contrac- 
tions of T ,Jll/ (qi,q2) with the four- vectors k^ 1 and Af to 
determine the coefficients Cj, i.e. we need to calculate the 
(— , _L)-components of T fJ " v (qi,q2) in modified light-cone 
coordinates. Since fc M has only a +-component, the eval- 
uation of the derivatives is now considerably simplified. 
Leaving the term proportional to pp' aside, we see that all 
derivatives which act on kj or kj' can be ignored. They 
would result in the replacement of p^ or p'^ by k^. Since 
k^ v = k^ v = k^*^ = and k 2 = kAi = 0, we do 
not need to consider those terms. The derivatives acting 
on kj or kj' are therefore only important to determine 
the term proportional to pp' . However, this is achieved 
more easily if the calculation presented in section llllBl is 
repeated with a scalar source term Jpp' in the exponent 
(see section HlIE[) . 

To calculate the preexponent of the polarization oper- 
ator, we must now insert the explicit expressions given 
in Eq. (f!M|) into the trace in Eq. (l4"2"j) . Since contracted 
with k^ or Af from each side, many terms of the trace 
vanish, e.g. the terms proportional to S 7 *", : 5 r2/il/ , 3 2/ip ?V 
Using the relations in appendix [D] we can show that 
Eq. (|42l) can be substituted by the following expression 



raV+pV 



yy 



+ <T [G 3 p$ v p' + G 3 p$ xP ' - 2Gt(p$ 2 xy p') - (pp')} 

- g 3 [($ v p'rp v - {%vpYp ,v +pn$ x p'Y - p"*($xpy] 

- Gi [p^idyPT +P'^vP) V + (wv + (S.pTpI 
+ G\ [(Z x pT(Z vP r + (W) M (3VP)"] 

- gi [(dypn^p'r + (ZyP'nzxpr] , (70) 

where &Jg = f(b)J/(%) = 3^ (ky)$ p " (kx) [here 
the replacement p^ — > (—i)dj andp' M — > (— i)d% is un- 
derstood if the trace is inserted in Eq. (|6Tj) . see Eq. (|5Tl) ]. 
Since the term proportional to pp' enters as g^ v , it mod- 
ifies only the diagonal coefficients C3 and C4. 

D. Evaluation of the derivatives 

Leaving the term proportional to pp' aside, we can ig- 
nore derivatives acting on kj and kj' as discussed above 
[this implies that the derivatives do not act on kp, kp' 
and ky, see Eq. (j54[) ]. The remaining non-trivial source- 
dependent part of the phase is given by [see Eq. 



s+t 



[tqai - sq 2 f + |(j + f) 2 + h{]+f)A ■ (71) 



The squared term contributes only if both derivatives act 
on it, which results in the replacement 



p a p'P 



n2oqq/3 



\ V U 3 °'J' f 2(s+t)9 



a0 



(72a) 



and the only non-zero contribution comes from the first 
line in Eq. (f70|) . If the derivatives act on the other source 
terms, we obtain the replacement 



p a p'? 



{-ifdftf, 

-> -^(^2 + £* a )(«£ - &*)• (72b) 

After these replacements are applied to Eq. ((70)1 and the 
sources are set to zero, we obtain (effectively) the follow- 
ing expression for Eq. (|70j) 

+ 2-^^ly^-PP'] 

+ w^i^y^Yi^iY -v 2 (z x qnz yq r] 



(73) 



[note that terms proportional to (3A) M , (3^Y can be 
omitted] . By changing the proper-time integrations from 
s, t to t, v [511 



ds dt f(s,t) 



+ 1 <.oo 

du / d,TTf(r, v) 
1 Jo 



(74) 



we see that the terms linear in v vanish. Those terms 
determine the coefficients cq — eg, which are therefore 
zero (as already anticipated from Furry 's theorem). 



E. Scalar term 

To determine the term proportional to pp' , we add 
the scalar source term ijpp' to the phase (instead 
01 WnJ 1 * + ip' u j '' 1 ) an d repeat the calculation presented 
in section UlI Bl The propagators are represented in the 
same way [see Eq. ([50]) ]. and we replace pp' by —i-£j- 
Then we take the integrals in dx + , dx- 1 , dy + , dy- 1 , dp'~ , 
dp ,J - , dp' + and dp + . Instead of Eq. (f52|) we obtain now 



(2tt)%- -x~- 



4st-J^ 



2{s+t+J)^ J 

< (2tt)<5[2( S + t + J)p- + (2< + J)q~}. (75) 

The remaining part of the phase (including the part from 
the propagator) can be written as 

iS' T = i[q^y~ - qt x ~ ~ P ± P ± J 

+ (-p- L p- L -m 2 )(s + t)-p- L A- L +Al. (76) 



It is now convenient to shift the proper-time integrations 



8-hj, t 



t 



-7 
2 u ■ 



(77) 



Due to this shift also the integral boundaries of the 
proper-time integrations depend on the source. However, 
if the derivative acts on the integral boundaries, either s 
or t is set to zero or to infinity. The resulting terms do 
not dependent on the external field since s = or t = 
implies /i = 0, ky — kx and thus A^ = and A = 0. 
On the other hand the terms at s —¥ oo or t — > oo do 
not contribute because the field-dependent part of the 
integral is convergent. Since we want to calculate only 
the field-dependent part of the polarization operator [see 
Eq. (f89|) ] we will ignore the source-dependence of the in- 
tegral boundaries. 

After the shift in Eq. ([77)1 . the delta functions read 



(2tt)%- 



(2n-J)q-]x 

x (2n)6[2p-(s + t) + 2q-t\ (78) 



and the phase is given by 

iS' T = i [(q£ - qt)x~ + (/i - \j)ql -m 2 (s + t- J) 
-p x p x (s + t)-p ± ^ ± +A\. (79) 

We can now use the delta functions to take the integrals 
in dy~ and dp~ (we keep writing y~ and jT for conve- 
nience). We then obtain the identities 



kp = --^kq, 

ky — kx + 



kp' = ^- t kq, 
(2\x-J)kq 



(80) 



[for J = this agrees with Eq. ([M])]. The shift in the 
proper-time integrals has the advantage that kp and kp' 
are now independent of J. We could have proceeded 
similarly also in the calculation of the other terms. How- 
ever, since we ignored sources contracted with k anyway, 
this was not necessary. 

Taking now the Gaussian integrals in dp 1 , dp 11 , we ob- 
tain the prefactor ., .*. , and the final phase is given by 

iS' T = i[{q+ - q+)x- + (ji - \j)q\ - m 2 {r - J) 

+ T m 2 J2£( I ?- J i)]> ( 81 ) 

i=l,2 

where Ii and Ji are defined in Eq. (|M|) [for zero sources 
Eq. (I5T|) agrees with Eq. ([55])]. Since pp' in the preex- 
ponent is only multiplied by g^" [see Eq. ([75])]. the eval- 
uation of the derivative is not very cumbersome and we 
obtain the replacement 



pp 



H) 



a 



\<& 



+ ™ 2 iE & Wi (ky) - 2 JiV* (ky)] (82) 



(84) 



after J is set to zero (as explained above, we have ig- 
nored the source-dependence of the proper-time integral 
boundaries). 

To symmetrize the final expression, we change the x~ - 
intcgration by defining [see Eq. (|62p ] 

z- = x~ + (pt - \J) q~ (83) 

(z~ = z~ for J = 0), which implies 

kx = kz — (fi — ^J)kq, 

ky = kz + (fj, — \j)kq 
and 

(q+ - q+)x~ + (M - \J)<& 

= (4 - lt)z- + (/! - \j)qiq2- (85) 
Finally, we obtain the symmetric replacement 

pp' — > (-i)^j — > -\q\qi +m 2 + m 2 ^ J^ & x 

»=1,2 

x [|V?(fcc) + H ? (%) - hMkx) - hA(ky)] (86) 

(we assume that at x~ = ±oo the external field vanishes 
and therefore the derivative does not act on the integral 
boundaries, which now also depend on the source). 

F. Final result 

To determine the non-vanishing coefficients c\ — C5 of 
the polarization operator [see Eq. (|68[) ] we combine now 
Eqs. dHJ), flSU), 03) and ([SB]). Furthermore, we define 
the following functions 

X tj = [I i -il>i{ky)][J i -il> j (kx)], 

z t = \[i)i(kx) - i>i(ky)] 2 

and note that for j = f = 

$%»A iu = -fk^^(kx), 
A" = -2mr £ Af&Ii, 

i=l,2 
eAi^^q^ = mkq^ipi(kx), 



A 1 A = 2mT&I t , 
eq$ x "k = 2kqTm 2 ^ t%ipi(kx)Ii 



i=l,2 



e 2 q$l y q = m 2 (kq) 2 £ tf^hx^ky). 



Using these relations, we obtain for the field-dependent 
part of the tensor T Ml/ the following expression 

T^( gi ,q 2 ) - T£: (g 1)<72 ) = ^vre 2 ^--^(fc - <z 2 ) X 
x do — dz- [6iA^A^ + 6 2 A^A? 

J-l ./o r J-00 

+ fe 3 A? A? + 64 A£ A^ + 65 Q£ QS] e l * , (89) 
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where the field-independent phase reads [see Eqs. (J59I) 
and |(B3J|] 

e J * = exp {i [{q+ - q+)z~ + mi q 2 - rm 2 ] } (90) 
\fi = \t{1 - v 2 ), see Eq. ||5B]>] and 



^T A 2i) e 



4^' 

b 2 ^2m 2 ^ 2 (^X 21 -^X y /), 
b 3 = - (i + *&) {e iT ? - 1) 



ir/3 



2m z 



±(eiZi+e2z 2 )+eiXn 



(i + ^)( e ^-i) 



2m^ 



^(^Z 1+ e 2 2 ^ 2 )+^22 



Jt/3 



p »T/9 



65 = -f (e JT " - 1) • 
The field-dependent phase is given by [see Eq. (|59[) ] 



(91) 



e^ 3 = exp [*Ym 2 ££(/*- J,)], 

where [see Eq. (64)] 



(92) 



8=1,2 



i; = 5 / d\ipi(kz — X/jkq), 



J, 



and [see Eq. {57}] 



r+l 
i / d\ipf(kz — X/ikq) 



(93) 



X; 



[I, - ^(fcz + /xfcg)] [Jj - ^(fcz - nkq)], 



Zi = h[ipi(kz - fJ-kq) - ipi(kz + fikq)Y 



(94) 



We note that, using the metric tensor g^", we can con- 
struct the following projection tensor [52j 



G^ 



q 2 ii -qmg 1 ^ , 



which obeys 

q^G^ = G^q 2v = 
and can be decomposed as 

GT = <zig 2 (A? A? + A£A£) + Q^Q". 



(95) 
(96) 
(97) 



This shows that the decomposition given in Eq. (|M|) has 
the structure claimed in 1521]. 



IV. DISCUSSION OF THE RESULTS 

A. Comparison with the literature 

The expression we obtained for the field-dependent 
part of T^ v in Eq. (|89p is manifestly symmetric in q\ 



and q 2 . We will now show how the alternative represen- 
tation found in [51] can be derived from our calculation. 
To this end we do not apply the shift in Eqs. (|62[) and 
(j83[) which means that we have to use the replacement 
given in Eq. (J82I) [rather then Eq. (|86|) ] for the pp'-tevm 
in Eq. (f73|) . which modifies the coefficients &3 and b^. 
Furthermore, we introduce the variable 

z'~ = x~ + 2[iq~ = z~ + /J<7~, (98) 

which allows us to write [see Eq. (|63|) ] 

($2 - lt) x ~ + wl = (it - <li) z ~ + Wi<72 

= (q+ - q+)z'- + fiq 2 (99) 

and [see Eq. (jM)) ] 

kx = kz' - 2fikq, ky = kz' (100) 

(here and in the remaining subsection we assume that all 
sources are set to zero). Thus, we obtain the following 
representation [see Eq. (J6TJ]) ] 



Ii = / dX ipi(kz' — 2kqfiX), 



I 2 - J 



i] 



Ji= dXip 2 (kz' -2kqfiX), 
Jo 

1 ,o /.l 



(101) 



dXA^/iX) 



dXA 2 (^X), (102) 



where we introduced 51 



Aj(r) = 4>i{kz' - 2kqr) - V* (&*')• ( 103 ) 

Furthermore, it is useful to define [compare with Eq. (|87l) ] 
AT, 



w., - [ii - V>;(%)] [I j ~ iJ>j(kx)], 
Yi = [Ii - ipi(ky)} [ipi(kx) - ipi(ky)] 



(104) 



which can be written as 

»i 



Xij - 



Yi 



o 



dXAii/iX) 



dXAi(/jX) 



dXAj(nX) - Aj(n) 



Ai(/i). 



(105) 



Finally, we obtain the following alternative represen- 
tation for the field-dependent part of T^ v 

T^( gi ,q 2 ) - T^ ( qil q 2 ) = -nre 2 ^-'^ - g 2 ) x 

r+ l r°° d T r+°° 
x / dv — dz'- [biA^AY, + 6 2 A£A^ 

J-l JO T J-OO 

+ b' 3 AlA u 1 +b' 4 AZAZ + b 5 Q?Q 2 J ]e i * (106) 
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with the coefficients 



61 = 2m a 66 (j-,X 12 - i£x 21 ) e*?, 



it/3 



b 2 = 2m 2 66 (^21 - ^A 12 

6i = - (7 + 4) (^ ~ 1) 



2m' 



(107) 



»i = - U + f ( 



pJT/3 



4 M 

-V (jT/3 _ 1 



(gYi + £y 2 ) + £x 22 



+ 2m z 
h = -& (e* 
and phases 

e<* = exp {« [(q+ - q+)z'- + m { - rm 2 ] } , 
e*' = exp[mn 2 ££(#-■*)] ■ 



(108) 



4=1,2 



This representation coincides with Eq. 2.27 in [5l|. 



B. Constant-crossed field 

The polarization operator for a constant-crossed field 
was first obtained in (z3, Iz3l (see also ® iHlzl). We 
show now how this result can be obtained from the ex- 
pression in Eq. (|89|) . 

A constant-crossed field is characterized by 



M4) = <t>, M<f>) = 0, 



(109) 



(the latter condition corresponds to £2 = and we will 
write £ = £1 in the following). The field tensor and its 
square are then given by [see Eq. (|11[) ] 

m 2 f 2 
e z 
For a constant-crossed field we obtain 
h = kz, J x = (kzf + \{pkq) 2 , I 2 = J 2 = 0, 

111 = -{ixkqf, Z x = 2(pkq) 2 , (111) 

Z 2 = X\ 2 = X 2 i = x 22 = 0. 

After inserting these expressions into Eq. (|89[) . we can 

take the integral in dz~ and obtain the remaining delta 
function 2w5^ + '(qi — q 2 ) 1 which implies that the polar- 
ization tensor for a constant-crossed field is diagonal in 
the external photon four-momenta. We define therefore 
the four-vectors [see Eq. (|67[) ] 



•jA 1 — nV — rt<" 



«? = <£, &> = (% = (% = 



k^q 2 — q^kq 
kq 



They obey 



kQ=-kq, qQ = 0, Q 2 = -q 2 . 



(112) 



(113) 



The four- vectors q^ , Q^ , A^ and A2 form a complete set 
and we obtain the following representation of the metric 
tensor 



1 



r 



{q^q v -Q' i Q v )-KlMi-K^K v 2 . (114) 



From Eq. (|89p we obtain now the following representation 
of the field-dependent part of T^ v in a constant-crossed 
field [see Eq. (ITTO)) ] 



T^( gi ,q 2 ) - T^ ( qi ,q 2 ) = -2iv 2 e 2 S 4 ( gi - q 2 ) x 

r+l 

x / dv I — [& 3 A^ + 6 4 A 2 1 A^ + 65Q^Q ,y ]e J * 



(115) 



where 



6a = -(i + 4)(e^-l) 

+ m \ 2 T 2 \{l-v 2 )[l-\{l-v 2 )]e^, 

h = -(7 + i) (e lTf3 - 1) + m 6 X 2 r 2 |(l - v 2 )e^ , 
fe 5 = -I(l~« 2 )(e^-l) 



and the phases are given by 



z$ = —ira, 



1-|(1-- 2 ) 



2x 9" 



*T/3=-f-r 3 &, NmY[l(l-» 2 )] : 



(116) 



(117) 



(in the following, we will make the change of variables 
r — > t, where r 3 b = t 3 and p = -§tt)- Here we have 
introduced the quantum non-linearity parameter 



X = 



eJWH _ ,V(kqr 

3 S " 



m - 



(118) 



We can rewrite now 



{FqYJFqY 
(Fq) 2 ' 
(F*qY{F*qY 



A 2 A 2 - {F * q)2 , 

where 

(F*q) 2 = (Fq) 2 = -^(kq) 2 

and obtain [see Eq. (j95|) ] 

G m^ = ^ _ g 2 ^ 

^(A^ + A^ + Q^ 

We note the following relations 

q p G pv = G""?, = 0, 
^G"" = G^kp = -fcgQ", 



(119) 



(120) 



(121) 



(122) 
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G W F 2 G av = ^ 2 {kq) 2 QYQ v . (123) where 



To obtain the representation given in |48l 15GJ, we pass 
over to different basis tensors 



fe 3 A^ + 6 4 A£A£ + b^Q!> = (q 2 b 5 - b, 



,2 U uAF*qy(F*qf 



+ (^65-64)- 



(F*q) 2 
and define the following functions [48|, [5C 

poo 

f(x) = i / dtexp [~i(tx + ^t 3 )] 



(Fq)»(Fq)» 

iJ (Fq) 2 

b 5 G" u (124) 



= 7rGi(x)+i7rAi(x), (125) 

pOO 

f'{x) = tdtexp[-i{tx+±t 3 )], (126) 

Jo 

f°° dt 
h{x) = / — exp (-*te) [exp (-ft 3 ) - 1] 
Jo l 



dt 



/(*)" 



(127) 



and 



/■°° ^ 
Ms) = / ^ exp (-ite) [exp (-f f 5 ) - 1] 

= -i[xf 1 (x)+f'{x)], (128) 

where Ai and Gi denote the Airy and Scorer function, 
respectively [76[ • These functions obey the following dif- 
ferential equations 



/"(as) = xf(x) - 1, 

/ x (x) = i-/(x) = -£/"(:*). 



(129) 



Using the latter and assuming convergence and vanishing 
of boundary terms, we can replace the function /i(x) by 
f'(x) in the following way 



+ 1 r+l 

dvg(v)fi[p(v)} = dv 

-l J-i 



G(v) 



p{v) 



f[p{v% 

(130) 



where G'(v) = g(v). 

Using the above notation, we can represent the field- 
dependent part of the tensor T M " for a constant-crossed 
field given in Eq. (|115|) by 



T^(quq 2 ) - T^ ( qi ,q 2 ) = i(2n)^(q 1 - q 2 ) X 

{FqY(FqY {F*qY{F*qY n 3 y 

n l /rn X9 1" ^2 — ~ ~G^ 



(Fqf 



{F*qf 



(131) 



e 2 m 2 r +1 



V3 



"= !^/_T* <TO+2) (»f m (i32) 

n= .pir dv m 

Air tt J_ 1 w 

[& = l(l ~ « 2 ) ; P = (f ) 2/3 (l ~ ££)]• Since all non- 
vanishing functions are even in v, we can now apply the 
following change of variables 



1 r 1 r°° 4 

dv = 2 / dv ~ dw — , 

i Jo J 4 w^w(w — 4) 



(133) 



which shows that the result in Eq. (|131[) is equivalent to 
the one given in [48|, [5(| . 



C. Quasi-classical limit 

We consider now a linearly polarized plane- wave field 

Vi(0) = V#), M<f>) =0 (134) 

(£ = £i, fi w = f^ u ) in the quasi-classical limit defined 
by £ ->• oo while [see Eq. (fTT8|) ] 



eVqPq _, V{kq^ 

%■ 3 ? 2 



(135) 



is kept constant. In the optical regime (photon energy 
cj ~ 1 eV) x > 1 requires £ 3> 1 which means that 
the quasi-classical limit is in general sufficient to analyze 
strong-field experiments with optical lasers. 

By employing the identity \kq\ = m 2 x/£,, we can ex- 
pand all functions depending on [ikq 

ll-Jl = -(l/i){pkq) 2 [^{kz)\ 2 + 0{pkqf, 

Z 1 =2(pkq) 2 [ij'(kz)] 2 + 0{fikq) 3 , (136) 

111 = -( f lkq) 2 [^'(kz)] 2 + 0(flkq) 3 

(X 12 = X 21 = X 22 = Z 2 = h = h = for lin- 
ear polarization). Thus, if multiplied by £ 2 only the 
leading-order terms are independent of £ and all others 
are suppressed. In the limit £ — > 00 the expressions in 
Eq. (J136I) correspond to those in Eq. (Illl[) with the re- 
placement x ~~* x(kz) = x^P'ikz). The remaining cal- 
culation is therefore similar to the constant-crossed field 
case and the final result in Eq. (1138)) corresponds essen- 
tially to Eq. (I131[) with the above replacement. Using 
[see Eq. ([TTTJ)) ] 



A^A* 



A£A£ 



(fq) 2 ' 

(f*q) 2 



(137) 
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and Eq. (|97p . we obtain for a linearly polarized plane- 
wave field in the quasi-classical approximation the fol- 
lowing representation for the field-dependent part of the 
tensor T^ v [see Eq. 



T^(q h q 2 ) - T^ ( qil q 2 ) = z(2n) 4 6^ ± \q 1 - q 2 ) x 



i /*+oo 

x — / dz~ e i{ ^- q ^ )z ~ 

2tt 



, U1YU1Y 



1 (fi) 2 
, {tiYiriY A 



where [see Eq. (fl"3"2"j) ] 



if*q) 



— —— / dw(tw - 1) 

47T 37T ' 



<7i 92 



■G^ 



, (138) 



e 2 m 2 ' +J 



lx(**)l 



-,2/3 



9 2 

e m 

47T 37T 



e gi92 

47T 7T 



dv (w + 2) 



-l 



w 

|x(Ml v: 



/'(p), 

/'(p), ( 139 ) 



! W 



with £ = K 1 - y2 )' p = .--7- 

G"" = q%q\ - qiq 2 g^ [see Eq. (gSD 



2/3 






Thus, we can write the field-dependent part of the 
tensor T^ v for a circularly polarized plane wave as [see 
Eq. (S3)] 



T^(q u q 2 )-T^ ( qi ,q 2 ) = -*7r e 2 ^--^fa - g 2 ) x 

x dv — dz- [b+A^A^ 

J-l JO T J-oo 



+ b-AtA v _ + i(6 x - & 2 )(A^ - A£A?) (W6) 

+ I (6 3 + M(A?A? + AgAJ) + fosQ^Q^e 1 *, 

where we defined 

A^ = A^ ± iA% (147) 

and the coefficients are given by 

b± = \ {(b 3 - b A ) T »(6i + b 2 )} = |m 2 e 2 x 
x [sinc 2 (Aifcg)-2sinc(2^fcg) + l] e T2ikz+iTP , (148) 



I(6 1 -6 2 )=m 2 C 2 {i±g|[- sm(2f,kq) 

+ 2 sinc(^fcq) sin(^fcq)] e 4T ' 3 , (149) 



D. Circular polarization 



The general result in Eq. (189)) also simplifies con- 
siderably if the plane wave is circularly polarized and 
monochromatic 

lM$ = Re f '*, ^2^) = ^*, 6=6=^ (140) 
We then obtain 

/1 = sinc(^fcg)3?e lfez , J 2 = sinc(^fc<j)3e lfez , 
Ji + J2 = 1, Zi + Z 2 = 2 sin 2 (/zfcg), 



(141) 



(142) 



(143) 



h ~ ipi (kz + fikq) = -ft A, 
h - ip 2 (kz + fikq) = $$A, 
h ~ ^\(kz — fikq) — 5R_B, 
h - ip 2 (kz - fikq) = S-B, 

where 

A = e lkz [sinc(fikq) — cos(fikq) — i sm(fikq)} , 
B = e lkz [sinc(fikq) — cos(fikq) + i sva(fikq)] . 

[we define sine a; = (sinx)/x]. Thus, 

X 12 -X 21 =SA*B, X n -X 22 = ftAB, 
X 12 + X 21 = SAB, X n + X 22 = ^A*B, 

where 

A* B = sine 2 (fikq) + cos(2 fikq) — 2sinc(2fikq) 

+ i [— sin(2fikq) + 2 sinc(fikq) sm(fikq)] , (145) 
AB = e 2ikz [sine 2 (fikq) - 2 sinc(2 fikq) + l] . 



l 2 (h + b i ) = -(± + ^)(e"f i -l) 

+ m 2^r 2 (i+V) sin 2 ( k } 

S l (i-v-) U* V (15Q) 

+ sine 2 (fikq) - 2 sinc(2 fikq) + l]e jr/3 , 



h = -^ (e"' 3 1) 

and the phases read 

ir/3 = irm £ [sine (fikq) — l] , 
z$ = i [(q% - q~l)z~ + miq-2 - Tin 2 ] 



(151) 



(152) 



[fi = \t(1 — v 2 )]. Finally, the integral in dz~ can be 
taken and we obtain the following expression for the 
field-dependent part oiT^ v (qi,q 2 ) for a monochromatic, 
circularly-polarized plane- wave laser field 

T^(qi,q 2 )~T^ ( qi ,q 2 ) = ^2~ dv 

/•OO 1 

/ - [T ^,5 4 ( gi - q 2 ) + T^5\ qi -q 2 + 2k) 
Jo T 

+ T^S i (q l -q 2 -2k)]e^, (153) 



l i44 i where 



i$ cp = -irm 2 {l + i 2 [l - sine 2 (fikq)]} +ifiqiq 2 , (154) 

T[T = tx(A»A v 2 - A^) + r 2 (A^ + A^A 2 ^) 

+ T3Q1QI (155) 
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T ± = 3^ 2 [sine 2 (jikq) 



2 siac(2fikq) + l] A^ A^ (156) 



and 



Tl = m 2 e^9)[2 sin 2 (nkq)/(nkq) - sin^fcq)] , 
r 2 = 2m 2 e £±g sin 2 ( M fcg) 



T3 



+ [(£ - |)«i*» - m 2 ] (1 - e"^) 



(157) 



This result agrees with Eq. 2.34 in [5i(. The terms de- 
scribed by T± u can be interpreted as describing processes 
where two photons from the background field are ab- 
sorbed or emitted, respectively (since the external field 
is not quantized, this interpretation relies only on the 
momentum conserving delta function). 

In order to obtain Eq. (|153[) from Eq. (|146[) we have 
used the identity 



r°° dr 

/ — e 4 *m 2 £ 2 [sinc 2 (Aifcq) - 2sinc(2^fcg) + l]< 
Jo T 



,ir/3 



dr 



/o T 
which follows from 



-d-(J T _ 1^ — i-± e iT P 



e^ + ^-m 2 ]^-!), (158) 



*^(^ -1) -',,-■ 

= m 2 £ 2 [sinc 2 (^fc<j) - 2 sinc(2^£;g) + l]e iTp (159) 

via integration by parts. 



CONCLUSION 



is measured in Heaviside-Lorentz units (eo = 1). The 
electron mass and charge are denoted by m and e < 0, 
respectively. Thus, the fine-structure constant is given 
by a — e 2 /(47r) ps 1/137. In covariant expressions the 
space-time metric g^ v with signature (1,-1,— 1,-1) is 
used and <9 M = (d/dt, V) is the four-derivative. This 
implies d^x v — g^ u , where x^ — (t,x) denotes the po- 
sition four- vector. The unit tensor is denoted by 5^ = 
g w 9pv = diag(l, 1, 1, 1) (5£ = 4) and space-time indices 
(lowercase Greek letters) are raised and lowered using the 
metric a p = gp V a v (summation over all types of repeated 
indices is understood if they do not appear on both sides 
of an equation). Greek and Latin indices take the val- 
ues (0,1,2,3) and (1,2,3), respectively. Contractions of 
four-vectors are denoted by a^b^ = ab, scalar products of 
three- vectors by a 1 b l — ab. We denote the dual of a sec- 



ond rank tensor T» v by T* 



{IV 



l e ^pa T where ^upa 



is the totally anti-symmetric tensor in four dimensions 
with e 0123 = — eoi23 = 1- For contractions of second- 
rank tensors and vectors a matrix notation is sometimes 
used, e.g. aTb = a^T^b u , (T^)^ = T?T% V , T 2 ^ = 



rppprp 



■ , (Ta)^ = T^ v a v . All spinors are Dirac spinors 
(with four components), spinor indices are usual sup- 
pressed. The Dirac gamma matrices are denoted by 7 M , 



(7 M 7" + 7V 



-«7°7 1 7 2 7 3 



and 2a^ v 



ryPryt* 



2g tiP ). For a spinor u we define u — u'7 






and for a matrix in spinor space M correspondingly 
M = 7°M^7°. A quantization volume V = 1 is as- 
sumed for the normalization of the single-particle elec- 
tron, positron and photon states. The total derivative 
of a function with respect to its argument is denotes by 
a prime f'(x) = -4-f{x). Integrals without boundaries 
range from —00 to +00. We use iO as a short nota- 
tion for ie together with the limit lim e _ >0 +- I n general, 
our notation therefore follows [55j with different units for 
charge. 



In the present paper we have proven for the first 
time the Ward-Takahashi identity for general loop dia- 
grams in a plane- wave background field (see section ITTEl) . 
Moreover, we have presented a new derivation of the 
leading-order contribution to the polarization operator 
in a plane-wave background field for arbitrary polariza- 
tion and dependence on the plane-wave phase (see sec- 
tion |TlT|. Our calculation relies on a direct evaluation of 
the space-time integrals without using Schwinger's op- 
erator method [4| that was employed in [51] . Finally, 
we have also shown explicitly why many coefficients of 
the polarization operator vanish [see Eq. (j6"5|) ] [5 ill I52J . 
The interesting feature of our final representation is the 
manifest symmetry with respect to the external photon 
four-momenta q\ and q 2 [see Eq. (|89l) ]. 



Appendix A: Notation 

In this paper we use natural units h = c = 1 (in some 
formulas H and c are restored for clarity) and the charge 



Appendix B: Light-cone coordinates 

Calculations involving plane-wave background fields 
become particular transparent if light-cone coordinates 
are used [68|, UM UM- Since the non-trivial space-time 
dependence of the momentum-space vertex in Eq. (|24[) 
is due to the plane-wave phase <f> — kx, it is natural to 
work in a basis where fc^ is one of the basis four-vectors. 
However, since k 2 = 0, this will be a light-cone basis. We 
introduce now a general light-cone basis by adding three 
four-vectors k^, ef (i £ 1,2) to the set and require the 
following orthogonality relations 



kti = ket — 0, kk 



-L, &i&j 



and the orientation 



tpvpijk k E1&2 — !■ 



(Bla) 



(Bib) 



To be more specific, we can in a reference system where 
the plane wave propagates along the direction n take the 
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following four-vectors 



k»=w(l,n), k» = —(l,-n), 

2u> 



(0,e 4 



ei x e 2 



(B2) 



(ei represent the two polarization directions of the plane- 
wave field and u has the dimension of a frequency) . 

Due to the relations given in Eq. (|B1[) . we obtain the 
following decomposition of the metric 



9tiu 



k^k u 



k^ku 



^lfi^lu — 62 M e2i/- 



(B3) 



This allows us to define the transformation to light-cone 
coordinates (primed indices) by 



AV":V 



b,A- 



A 



-ip 



•A" =6', (B4) 



where the components denote the following scalar prod- 
ucts 



A M — &A" 



A V = - 



A" 



£• A 11 



e2 M 



(B5) 



Since all properties of the light-cone coordinates follow 
from the relations in Eq. (jBip , we are not forced to use 
the canonical basis in Eq. (|B2|) . For the calculation of 
the polarization operator it is more convenient to use the 
two four- vectors [see Eq. (|45|) ] 



?'/' = A? = 






JH 



kq^/~- 



K 



J 2 1v 



kq^f- 



(Blla) 



together with fc M and 



fc'^ = k^ + 






a29 -, 



(axq) 2 (a 2 q) 2 



2{kqf 



k^. (Bllb) 



The set of four- vectors fc M , k ,fJ -, e'/ 1 , e^ also obeys the 
relations in Eq. (|B1|) and we will call the coordinates, 
following from this set, modified light-cone coordinates 
[the same symbols (+, — , i, n) are used to denote the 
corresponding components] . 



Appendix C: Gamma matrix algebra 



(we label light-cone components by +,— ,i,ii). On the 
other hand the inverse transformation is given by 



A- 1M 4 
a -im 



= fc", A 
= jfe", A 



-ip 






(B6) 



We point out that fc M has dimension of momentum and 
therefore k^ must have dimension of inverse momentum 
(ef are dimensionless) . Hence, the dimensions of v + and 
v~ differ from those of v^ (here v 11 is an arbitrary Lorentz 
four- vector) . The different dimensions of the light-cone 
components can be circumvented by defining k^ = un^ 
and using the dimensionless quantity n M in place of k^. 
Then, however, nv is not a Lorentz scalar (in contrary to 
kv = v~) and u> has to appear explicitly in many places. 
In light-cone coordinates the metric is given by 



9ii'V 



^A-^A" 10 



^ + S~ t S+, 



8^81 



8 n ,8 11 , 



(B7) 



which allows us to write the scalar product of two four- 
vectors as 



a^ = a + b- 



a~b + -a l b l -a ll b l 



(B8) 



(we also use the short notation a^b 1 - = a I b l + a ll b 11 ). Due 
to Eq. (jBip we obtain 



detA'', 



= |A+ A~ A 1 A 11 ^ vpa \ = i 



(B9) 



Thus, the four-dimensional integration measure becomes 
/ d A a = I da+da-da^, da x = da 1 da 11 . (BIO) 



In this appendix we summarize some general identi- 
ties, which are useful in calculations involving gamma 
matrices. The gamma matrices form a complete set in 
the sense that any matrix in spinor space can be decom- 
posed according to [79| 



r = cil + C57 + c^ + C5 M i7 M 7 + c^ia 



uv 



(CI) 



where we can assume that c M „ = — c vfi and the coeffi- 
cients can be calculated using 



ci = \ tr ir, c 5 = \ tr 7 5 r, c M = \ tr 7M r, 
c 5fi = itr«7 AI 7 5 r, c llv = \tria llv Y. 



(C2) 



Due to the cyclic property of the trace one can recursively 
calculate traces of arbitrary length without conceptual 
difficulties by permuting the first gamma matrix to the 
last position. For completeness we note the following 
relations 



itryyw 7 

l 

4 






ltr 7 Y7YT 5 = ie^ p(T . 



(C3) 



Thus, any identity involving gamma matrices can be 
proven by calculating the fundamental terms given in 
Eq. (IC1|) for both sides of the equation. It is in par- 
ticular possible to map the gamma matrix algebra to a 
corresponding tensor algebra once the decomposition of 
the product of two (arbitrary) gamma matrix expressions 
is known 



r„i\ 



(C4) 
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Here T x is written as in Eq. (jCll) with the letter c replaced 
by the letter x appearing in the index. The coefficients 
of r c are then given by 

c ± = ai&i + a 5 b 5 + a p b^ + a p b 5fl + 2a fny b fiU , 

c 5 = (0165 + 0561) + (ia^fofj, - ia$p,b p ) 

Op = {a-ibfj, + a,J>i) + (ia 5fl b 5 - ia 5 b 5fJ ,) 

+ 2 (ia au b u - ia v b^) - ie aup(T (a»b p ° + 0?°%) , 

C5p. = (aib 5 p + a 5M 6i) + (ia 5 bp - ia^bs) 

+ %^ vp „ {.a v b p ° + a p °W) + 2 (ia uv b v 5 - ia^) , 



Cfiv = (aib^ + a IMU bi) 



2 \*in>po 



(a pa b 5 + a b b pa ) 



- 5 (a M &„ - a u bf,) - ^ vp o {a p b% + of 6 P ) 

- § {a^b bu - a bu br p) + 2i (a w b p w - a up b p ^) . 

(C5) 

We point out that taking the trace of the gamma ma- 
trix expression T c projects out the coefficient ci [see 
Eq. (|C2j> ] . Therefore, one can also use Eq. (|C5|) in the 
calculation of traces. 



Appendix D: Tensor relations 



If Eq. (|C5jl is used to simplify large gamma matrix ex- 
pressions, one typically encounters products or contrac- 



tions of the totally anti-symmetric tensor e a Pi s . They 
can be simplified using well-known identities stated here 
for completeness [66J 



e a ^ S e aM& = -24, 



^ €a{3pcr 



m 



-S^6 p a + S p 6;S p a ~S p S^ p ), 



(Dl) 



^otp-yS 



= det 



/ *2 % % % \ 

K % «; ^ 

*s ^ a? <*? 

V *S Si 5° 8% J 



In particular, we note the following formulas for anti- 
symmetric tensors T Ml/ , T^ v and T% 



rp*fJ,Urp*Otp 1 

J l J 2 — 2 



f^V" 



9 9 ) i \p<r 1 2 - J l J 2 



g»«( Tl T 2 )PH - g^^Tz) 



0V 



9 up (TiT 2 ) ap + gW{T x T 2 ) av 



pPi 



(T-jT*V'" 



„> . - y^n^rf 



(TxTa)"'*, 



1 lp,v 1 2 ~ ~ J lp^ 1 2 



and 



e 

2 e IJ.vpiT 



fiuptrrp* XPrpUp _ XHrppp 

act u a w a 



. S P T V 



1 T*pff _ _T^ 

1 tHl/0(7 •*■ -*■ / 



(D2) 
(D3) 

(D4) 
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